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A mathematical model was formulated to describe the performance of a 
hydrogen-bromine fuel cell. Porous electrode theory was applied to the carbon felt 
flow-by electrode and was coupled to theory describing the SPE system. Parametric 
studies using the numerical solution to this model have been performed to determine 
the effect of kinetic, mass transfer, and design parameters on the performance of the 
fuel cell. The results indicate that the cell performance is most sensitive to the 
transport properties of the SPE membrane. The model was also shown to be a useful 
tool for scale-up studies. 


INTRODUCTION 

Hydrogen-bromine fuel cells are of interest as both primary and regenerative 
energy storage systems. The fuel cell system can be coupled with solar cell arrays 
to provide the power necessary to charge the system. The regenerative capabi lity of 
hydrogen-bromine fuel cell systems, along with their high energy densities, make them 
excellent candidates for space power applications. 

The electrochemical reactions for the hydrogen-bromine system are nearly 
reversible, and the use of solid polymer electrolytes has eliminated cell gaps. 
Consequently, good energy storage efficiencies can be obtained even at high current 
density operation. 

High current density operation requires reactor designs which enhance mass 
transfer rates. A flow system is used to improve mass transfer and to aid in thermal 
management . 

A promising design for a hydrogen-bromine device is one having a negative 
half-cell with only a gas phase which is separated by a solid polymer ionic 
conducting membrane from a positive half-cell with a flowing aqueous electrolyte. 

The hydrogen and bromine are stored external to the cell, the hydrogen in the form of 
a metal hydride. This configuration warrants theoretical treatment for further 
development and scale-up. Therefore, this work has focused on the development of a 
theoretical model capable of predicting overall cell performance. 

In this paper, a theoretical model of a hydrogen-bromine fuel cell is described. 
The model equations were solved by numerical techniques. The parametric studies with 
this model elucidate the effects of catalyst activity, mass and ion transport, and 
cell scale-up on single cell performance. A discussion is presented on the results 
of these studies. 


*Th i s Research was supported by NASA Con trac t. NAG 3-500 and NSF award CPE-8351849 
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SYMBOLS 


A = characteristic area of an SPE consisting of exposed and unexposed membrane, cm 
a = specific electrode area (i.e. the active area per unit electrode volume), cm * 

g 

aJ, = surface bromine activity 
br 2 


f 5 
$r. 


= bromine activity of reactor feed 


= local bulk bromide concentration, mol/1 

Br- 


C Br- 


= surface bromide concentration, mol/1 


Cj^_ = bulk |3-phase concentration of protons in equilibrium with the 6-phase, mol/1 




= proton concentration in the p-phase at the electrode surface, mol/1 


.f6 


Cgr- = cell feed concentration of bromic acid, mol/1 


= diffusion coefficient of bromide ion in solution, mol/1 
Br 

D u ^ = effective proton diffusion coefficient in the membrane, mol/1 
H+ 


F = Faraday’s constant, 96, 487 A-s/equil 

f = fraction of SPE not in contact with the current collector 

k = mass transfer coefficient, cm/sec 
m 

I = total current from a characteristics area (A^), A 

2 

i = current density in the porous electrode based on the true electrode area, A/cm 

i = bromine reaction exchange current density evaluated for a bromine reference 

2 

electrode at cell inlet concentrations, A/cm 

2 

i* = hydrogen side current density, (based on membrane area) A/cm 

i* r = exchange current density (based on membrane area) evaluated at unit proton 
o,rel 

2 

activity and hydrogen pressure. A/cm 
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L,L' = length and width of SPE section not in contact with current collector, cm 
n = number of electrons transferred in reaction, equiv/mol 

p u a = hydrogen partial pressure, atm. 

n - 2 

2 

R g = effective membrane area resistivity, ohm-cm 

t = porous electrode thickness, cm. 
tp = membrane thickness, cm. 

t = transport number for protons in the membrane 
U = equilibrium cell potential, volts 

v y = plug flow solution velocity through the porous electrode, cm/sec 

a . a c = anodic and cathodic transfer coefficients, respectively, and usually 
reported as both being 1/2 

6^ = thickness of platinum catalyst layer on SPE, cm 
e = void fraction of bromine electrode 

e' = volume fraction of membrane material between platinum catalyst particles 
jc = effective electrolyte conductivity, mho/cm 


o 

P = effective resistivity of platinum catalyst, p of solid platinum, ohm-cm 
e 0 

<f>l <#>2 = bromine side metal and solution phase quasielectrostatic potentials versus a 
bromine reference electrode at feed concentrations, V 

<f >2 = hydrogen half cell quasielectrostatic potentials relative to a normal 
hydrogen electrode of the metal phase and the /3-phase respectively, V 

THEORY 

The description of a single hydrogen-bromine cell is shown in figure 1. The 
hydrogen electrode is treated as a planar electrode having a current carrying bus net- 
work. The bromine electrode is a porous carbon felt. 

The planar hydrogen electrode, (a-phase), is in direct contact with the solid 
polymer electrolyte denoted as the /3-phase. The kinetic expression for the hydrogen 
reaction can be written according to the nomenclature and analysis of Trainham and 
Newman [ref. 1] as 
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( 1 ) 


i 


o.ref 


• a 
P H 

1 atm 


1/2 r F ^ 



1 mol/1 




where the nomenclature is presented above. 

The partial pressure of hydrogen at the catalyst surface is assumed to be equal 
to that of the bulk. The kinetic expression also requires specification of the 
concentration of the protons at the surface of the platinum particles. The Nafion 
membrane is believed to be composed of clusters of hydrophilic sites in a hydrophobic 
matrix [ref. 2]. Protons are associated with the fixed functional sulfonate groups 
and with the free water phase. For our model we assumed that the protons involved 
with the surface reaction come from the free acid phase. Protons are then 
transported across the membrane through a diffusion/migration mechanism. The 
interfacial mass transfer can be expressed as: 

i’ = ~ nFD H±_ [eg - C£ /] x 10" 3 (2) 

The concentration of protons in the membrane at the side opposite the hydrogen 
electrode is assumed to be in equilibrium with the adjacent HBr solution. This 
concentration was estimated by using a curve fit of the data given by Yeo and Chin 
[ref. 3], which is expressed as: 

= -0.458 + 0.80067C^ r _ + 0 . 023843C^ r _ 2 (3) 


We have recently begun work to examine the role of the protons associated with the 
fixed ions on the kinetics and to examine the proton diffusion mechanisms in Nafion 
membranes . 

The metal phase potential, (J)^, is taken to be uniform. A conductive bipolar 

plate feeds the current to the platinum catalyst. Although current must flow through 
a thin platinum layer to the bipolar plate contact point, the potential drop is 
estimated to be small. A simple model which justifies this assumption is described 
in the Appendix. 

The bromine electrode performance is described by macroscopic homogeneous porous 
electrode theory [ref. 4,5]. The necessary descriptive equations are now summarized. 
The charge balance, caui be written as 



(4) 


with boundary conditions 
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X 


(5) 


= 0 , = o 

dx 

x = t, <f> 2 = ^ + U - i’R s (6) 

where R g is the membrane resistivity. After neglecting diffusion in the axial 
direction, the conservation of mass of the bromide ion is written as 


ar 5 

Sir- _ 6 

y fly e °Br- 


9x 


al 1 r>3 

F X 10 

nF 


CO 


with boundary conditions of 


x = 0, t 


dx 


= 0 


y = °’ Sr- = C 


( 8 ) 

(9) 


The surface concentration is related to the bulk concentration through the mass 
transfer coefficient, k^, by the following expression. 


1 - nFk m [<t- - Sr-] * 10 ~ 3 

The kinetics of this reaction can be expressed by 


( 10 ) 


1 = 1 


cr- 


o. ref 


r . 5-, 1/2 




exp|^£(^ - ^2 exp|_^E.(^ -<l> )| 

'-RT 1 ^ ^ RT 1 2 f 




( 11 ) 


In these calculations, the ratio of the bromine activities is taken to be unity and 
the carbon phase potential, is assumed to be uniform. 

At the membrane-porous electrode interface, the membrane current density, i’, is 
related to the porous electrode current density, i, by the following expression: 


d<f> 


dx 


= i ’ = -a X i dx 
x=t 0 


( 12 ) 


The equations of this model were combined, made dimensionless, linearized, and 
then solved using an implicit numerical technique. The dimensionless equations are 
listed in Table 1, with the corresponding dimensionless groups given in Table 2. 

The numerical technique used here was that developed by Newman [ref. 6]. The 
nonlinear equations (B) and (D) of Table 1 were linearized by a first order Taylor 
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series expansion. Then the derivatives were written in second order finite 
difference form. The boundary at the bromine electrode-membrane interface was taken 
into account by two methods. In the first method (Equation (E)), the gradient of the 
potential in the solution phase is equated to the membrane current density. At the 
electrode entrance, the distributions are calculated using the boundary conditions 
and Equations (B) through (E) . A step is made in the axial direction, and Equations 
(A) through (E) , along with the appropriate boundary conditions are solved using the 
Cr ank-N i cho Ison symmetric form. All the equations are solved implicitly with Cj - C g 

unknown variables. The iterative technique was carried out until convergence with 
the Taylor series expansions taken about the solution of the previous iteration. 

An alternate method of handling the bromine electrode-membrane interface 
involved using equation (E') instead of equation (E) . The integration was carried 
out using the trapezoidal method. Since C 4 and C g are only needed at the membrane 

position, these variable were used as dummy variables for storing the summation terms 
of the integration across the porous electrode. This eliminated extra variables, and 
allowed the integration to be performed implicitly. 

The results of the numerical calculations were checked by comparing the results 
at the boundary conditions and by verifying intermediate calculations. The step 
sizes across the porous electrode and along the axial direction were varied to be 
sure that they did not have a significant effect on the results. The convergence 
criterion was increased from a relative error of 5 E-6 to 5 E-8. This had no effect 
on the calculated results. 

The calculations using equation (E) took ~7 iterations at the first position and 
about ~3 iterations at subsequent axial positions for the base case. We devised the 
second method (use of equation (E')) because it was thought to be more stable since 
high current-densities and consequently steep gradients were expected with this 
system. However, the number of iterations to convergence and computational times 
were about the same. Even when comparing the two methods with high overpotentials, 
the second method required a few more iterations and slightly more computer time. 
Therefore, any advantages offered by the second method were not recognized and the 
method did not appear to justify the added programming complexities which were 
necessary to implement it. 


RESULTS AND DISCUSSION 

Model calculations were first done on a base case of parameters which represent 
reasonable values for a small single cell. The values for these parameters are 
listed in Table 3. The feed bromide concentration corresponds to approximately a 17% 
state of charge of an initial solution concentration of 48% HBr. The membrane 
resistivity was estimated from data for a 24% HBr solution [ref. 3]. The effective 
solution conductivity was calculated from the conductivity of the electrolyte by 
[ref. 1]. 


1.5 

K = K 6 
O 

where k = 0.8Q _1 cm _1 [9] and a was estimated by knowing the density of carbon felt 

o 

[10] and the density of graphite [7]. The dimensionless axial length for the base 
case was chosen to be Y = 0.0282656. This length corresponds to a cell length to 
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v ® 1 ° c * ty 1 ' a ^ 1 ° of 77.5 s (or a 15.5 cm long cell operating with a solution velocity 
of 0.2 cm/s). J 

The open circuit potential, OCV, of a prototype SPE hydrogen-bromine cell was 
measured experimentally. The details of this study were reported elsewhere I" ref 111 
and only a summary of the findings will be given here. In figure 2 the mean ' 
experimental values of the OCV and the corresponding 95% confidence interval are 
shown as a function of state of charge. The OCV’s for solutions with a charge 
capacity of 48 wt% HBr (i.e. initial uncharged solution) are lower than the OCV’s for 
a system with a charge capacity of 35 wt%. This indicates that a trade-off exists 
between the OCV and the charge capacity of the system. 

The experimentally measured OCV’s of this work were compared to the 
semi-empirical correlation reported by Yeo and Chin [ref. 3], Their correlation was 
derived from experimental data assuming a unit activity of bromine. The correlation 
predicts somewhat greater open circuit voltages than those found in our experimental 
work at the lower states of charge. At the highest state of charge, the OCV 
predicted by the Yeo and Chin correlation was somewhat less than that found 
experimentally, even though the solubility of bromine in the experimental solution 
was exceeded (i.e. the bromine activity was unity). The differences could be due to 
the differences in tribromide and pentabromide concentrations. Overall, though, the 
comparison between our experimental data and the correlation of Yeo and Chin Tref. 31 
was good; See figure 3. Therefore, in this simulation study the correlation of Yeo 
and Chin [ref. 3] was employed for estimating the OCV of a hydrogen-bromine fuel 


Base Case Calculations 

Some calculation results for base case conditions are given in figures 4 and 5. 
The calculations were done for a charging mode of operation with a 200 mV applied 
polarization. Figure 4 shows the distribution of local current density across the 
bromine electrode. Under these conditions, only half of the electrode appears to be 
active. A larger current density occurs near the cell entrance. In figure 5 the 
distribution of bromide ion across the bromine electrode is shown. At the cell 
entrance, the bromide concentration is constant as established by the boundary 
condition. However, a significant concentration variation across the bromine 
electrode is noted when moving up the cell in the axial direction. Also, one can 
even note the drop in bromide concentration near the current collector. This, of 
course, is due to a net consumption of bromide ions. The significance of these 
calculations will become clear when variations on the base case are discussed. 

Bromide Ion Mass Transport 

The mass transfer coefficient associated with the bromine electrode for the base case 
was calculated from the Bird correlation (12) using a velocity of 1 cm/s, and 
assuming a cylindrical geometry for the fibers of the porous carbon electrode. The 
calculations show that there is very little effect on cell performance if the 
magnitude of the product ak m is lowered by one order of magnitude or increased by two 

orders of magnitude. Mass transfer of the bromide species to the carbon fibers does 
not appear to be a problem in this system. 
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Kinetic Effects 


The model calculations are useful for assessing the effectiveness of electrode 
catalyst. For the base case solution composition , the exchange current density was 
estimated [ref. 13] using data given by Mastragastino and Gramellini [ref. 14] for a 
smooth vitreous carbon electrode. As shown in figure 6 , an increase in the exchange 
current density will only have nominal impact on cell performance. For the hydrogen 
electrode, the base case exchange current density was estimated by assuming that the 
ratio of true platinum area to membrane area is 300:1 (i.e. an exchange current 

density of 0.3 A/cm ) . As shown in figure 7, a further increase in platinum loading 
or enhancement of surface area will not significantly improve cell performance. 
However, a loss in platinum activity could severely compromise the performance of the 
cell. This could occur, for example, if bromine species are transported across the 
membrane so that they adsorb on the platinum and thus poison the active area. 

Membrane Effects 

The model calculations demonstrate that the membrane transport properties have 
the most dramatic impact on cell performance. To demonstrate this, charge and 
discharge polarization curves are shown in figure 8. The performance is enhanced 
significantly when the membrane resistivity is lowered by an order of magnitude 
(psi/10). In this figure, a limiting current density is approached when in the 
charging mode. This effect is caused by a limit in the diffusion rate of protons 
across the membrane. The diffusion coefficient of protons was calculated by using 
the Nernst-Einstein relation and membrane conductivity data. If transport is 
restricted to only migration of protons associated with fixed ions, then a 
concentration gradient would not exist. However, since there is a significant amount 
of free acid in the membrane under these conditions [ref. 3], the transport number 
for proton does not necessarily equal unity. Thus, a concentration gradient can be 
established. Actually, the transport mechanism is likely to be more complicated 
since protons are moved by simple diffusion and migration in the free acid and by 
migration only along the fixed ion sites. Therefore, the effective transport number 
is probably not constant with current density. These results indicate that further 
detailed studies of membrane transport under these conditions are necessary. 

Another aspect of the membrane effect deals with the local acid concentration in 
the solution phase adjacent to the membrane being much lower than the inlet acid 
concentration. This was shown earlier in figure 5. The computer calculations were 
modified to account for this effect. The result is shown in figure 9 where the 
membrane current density along the axial direction is reported. Under these 
conditions, the changing acid concentration is shown to have a large effect on cell 
performance. 


Scale-up Effects 


Some scale-up considerations were examined using the model developed here. For 
example, increasing the length of the cell in the axial direction results in a lower 
overall cell current. This is shown in figure 9. Another question which often 
surfaces when designing a cell deals with ascertaining the optimum thickness of a 
porous electrode. The effect on cell performance of electrode thickness when the 
electrolyte velocity is maintained constant is shown in figure 10. As can be seen in 
this figure, under these conditions an electrode thickness greater than ~1.8 mm only 
provides marginal improvement in cell performance. As the electrode thickness is 
increased though, a large volumetric flowrate of electrolyte passes through the cell 
and consequently, the conversion per pass decreases. 
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On the other hand, if the volumetric flowrate is maintained constant when 
examining the effect of electrode thickness, the current density decreases as the 
electrode thickness increases. This is shown in figure 11. In this case, as the 
electrode thickness increases, the linear velocity decreases which also decreases the 
convective transport of mass. As a consequence of this, the concentration of bromic 
acid drop considerably which causes slower kinetics and thus lower current density. 

CONCLUDING REMARKS 

A mathematical model was formulated for evaluating design and operating 
parameters of a hydrogen-bromine fuel cell. The model coupled porous electrode 
theory applied to the bromine electrode with membrane transport and a planar 
treatment of the hydrogen electrode. 

The parametric studies demonstrate that both the bromine and hydrogen electrode 
specific surface area and kinetic parameters influence cell performance. However, 
the magnitude of performance enhancement decreases as catalytic and surface area 
effects increase. Mass transfer does not seem to be a significant limitation in the 
bromine electrode. 

The calculations indicate that the transport properties of the membrane are the 
most influential factor in controlling cell performance. Further research appears to 
be necessary to characterize the proton transport mechanisms in the membrane 
material . 

Finally, the model developed here was demonstrated to be a useful tool for 
studying cell scale-up effects. 
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APPENDIX 


Justification of the Equipotential Hydrogen Electrode Assumption 


An enlarged view of the porous platinum electrode with a current collector is 
shown in figure A1 . The electronic current must be conducted through the platinum 
layer to the current collector. The closer the current collector contacts are to 
each other, the lower the voltage drop will be within the electrode. However, many 
contact points will block the platinum surface and thus will cause higher 
overpotentials. This trade-off will now be examined. 

Assuming an approximate uniform current distribution, a voltage balance gives 
[ref. 15] 



LL' K e 86 


where L' is the width of the electrode section, and 6 is the thickness of the 

e 

platinum layer. In this expression, A(|)j is the potential drop in the platinum layer 

over the distance L/2. If f is the fraction of electrode area exposed to the gas 
phase, then for a symmetrical current collector system 

A t f = LL’ (A-2) 


where A^ is the total area of the membrane. The apparent current density can be 
written as 


i 1 



(A-3) 


Therefore, equation (A-l) becomes 
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(A-4) 


• p e L 2 


A( h- S L1L 
n f 6 8 


Of course, as f approaches unity, the contact area between the current collector and 
the platinum electrode diminishes and the contact resistance may become substantial. 
This effect, however, is negliected in this analysis. 

Assuming the platinum layer consists of a close packing of spherical particles, 
the Bruggeman equation can be used to relate the effective metal phase resistivity, 

p , to the resistivity of solid platinum, p °: 
e e 

Of. ,..-3/2 

P e = P e (1 ~ e > 

In this expression, e* is the void volume fraction between the platinum particles. 

A summary of calculations of the potential drop is shown below. These results 
indicate that the potential drop along the platinum electrode is small. Therefore, 
the assumption in the model of an equipotential electrode is justified in most 
instances . 


Calculation of AC^ Summary 



e' = 0.3 

p = 1.81 

e 

-5 

x 10 ohm-cm 


8 = 30 urn 

e 

r = iooo 

mA/cm^ 

f 

L, mm 

6 , pm 

A<j>*, mV 

0.2 

2 

5 

0.9 



10 

0.45 


5 

5 

5.65 



10 

2.82 

0.5 

2 

5 

0.36 



10 

0.18 


5 

5 

2.26 



10 

1.13 
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Table 1. Dimensionless Model and Boundary Conditions. 


t 

° 

II 

1 

1 

JO 

_C> 

i 

O 

S3 

(A) 

aY ax 2 


pc f (C 1 - C 2 ) = C 2 exp{(f) - C 3 } - exp{C 3 - <|>} 

(B) 

2 


^2 = -C f (C 1 - C 2 ) 

(C) 

2 

dJT 


1 /9 

"MCg - C 4 ) = P A/ exp{o - C 5 } - C 4 exp{C 5 - a) 

(D) 

dC 


^ C 4 - C 6> = d)T 

(E) 

1 


C 6 - C 4 = * C f f ‘ < C 1 - S)® 

(E ) 


Boundary conditions: 


o 

II 

(F) 

dX 


o 

II 

si 

(G) 

dX 

✓—s 

1 

u 10 

■9- 

+ 

CD 

+ 

II 

cf 

(H) 

O 

II 

si 

(I) 

dX 


C l = 1 

(J) 
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Table 2. Definition of Dimensionless Variables and Parameters. 


Variables 


c 1 = Sk,. 

c£ 5 

Sr- 


Cl 6 

C 9 = . 

P f 6 
C Br- 


C 3 = ff*2 

RT 




1 mol/1 


a . 

RT 




1 mol/1 


X = x/t 


Y = &D Br 
v t 2 

y 


Parajneters 


<t > = f££ . 

KT 


0 = gUF , 
RT 


or = _Ii_ , 
RT 


o.ref 


r f 5 

C, = Srz. 


1 mol/1 


p _ nFk (0.001 mol/cc) 


or = 


_ an F 2 at 2 k^(0.001 mol/cc) 


kRT 
Pi 


, nFD„7 (0.001 mol/cc) 

A = £1+ , 

t„i’ ,(l-t ) 

P o.ref + 


n = 


ale t 


eD, 


Br- 


a 

p = Ph 2 
1 atm 


= onF 2 D^R s (0.001 mol/cc) 
t RT (l-t + ) 


p 


D, 


H+ 


(0.001 mol/cc)aF 2 ntD^+ 


V = 


kRT (1 - t ) 
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Table 3. Base case values used in parametric studies. 

p^ a = 10 atm. ^Br- = 7-45 mole/liter t = 0.3175 cm. 

tp = 0.025 cm. [ref. 3] D R ^ = 9.49xl0 _7 [ref. 3] = 3.87xl0 _5 [ref. 7a] 

R g = 0.255 n-cm 2 [ref. 3] k = 0.74 Q 1 cm 1 [ref. 8, 9] a = 280 cm -1 [ref. 8] 

p = 1.49 g/cc p = 1 g/cm-sec k = 0.0866 cm/sec, 

m 

i o ref = 39.7 mA/cm 2 i’ Q ref = 300 mA/cm 2 U = 0.816 V [ref. 3] 

<*>. = 1.016 V e = 0.95 t = 0.975 
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Experimental OCV (volt*) 


Figure 3: Comparison of experimental open circuit voltages for solutions with 48% 

HBr charge capacity to open circuit voltages predicted by correlation of Yeo and 
Chin. 
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Figure 4: 


Local current density distribution in porous bromine electrode. 
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Figure 11: The effect of electrode thickness of cell performance with constant 

solution volumetric flow rate. 
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Figure A1 : Expanded view with dimensions of a hydrogen electrode of a hydrogen 

bromine SPE cell. 
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